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Abstract 

The  effect  of  the  excitation  conditions  of  a  traveling 
wave  tube  on  its  operation  is  analyzed  in  detail.  By  taking  the 
Laplace  transform  of  the  field  equations  with  respect  to  the  longi- 
tudinal coordinate,  the  initial  conditions  are  broti^ht  directly  into 
the  eq''JStions,  so  there  is  no  problem  of  matching  at  the  input.   The 
transformed  equations  iray  be  solved,  and  the  actual  fields  and  currents 
then  determined  by  inverting  the  transform.  The  current  in  the  ampli- 
fied mode  is  found  and  it  is  expressed  in  terms  of  transfer  coeffici- 
ents which  give   the  current  produced  down  the  tube  by  arbitrary 
distributions  of  input  field,  input  current,  or  input  a.c.  velocity. 
In  this  paper,  only  the  analytical  work  is  presented.  Applications 
to  specific  excitation  systems  will  appear  in  a  second  paper. 
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I.  Introduction 

1  2 
Previous  research  '  on  the  traveling  wave  tube  have  proceeded 

along  the  following  lines.  First  the  helical  portion  of  the  tuhe  is 
assumed  to  he  infinitely  long.  The  physical  helix  is  replaced  by  a 
cylindrical  sheath  which  is  conducting  only  in  the  helical  direction. 
The  properties  of  the  three  discrete  modes  of  this  system  are  then 
studied  in  detail.  To  find  how  these  modes  are  excited  by  the  input 
fields  of  the  act-'jal  physical  system;   the  input  electric  field,  elec- 
tron velocity,  and  a.c.  charge  density  are  then  expressed  in  terms  of 
the  coefficients  of  the  three  modes,  and  the  resulting  equations  solved 
for  these  coefficients.  This  theory  has  been  extended  phenomenologically 
by  introducing  a  parameter  into  the  equations  to  take  account  of  passive 
modes' ,  but  again  the  fields  of  the  finite  length  helix  are  replaced  by 
those  of  the  infinitely  long  system.  In  this  paper,  we  shall  study  the 
validity  of  these  processes  and  present  a  new  method  of  treating  the 
traveling  wave  tube  which  suiBmounts  the  difficulties  of  the  previous 
analysis. 

The  first  difficulty  faced  by  the  earlier  analysis  is  that 
the  discrete  modes  of  the  helical  sheath  system  are  not  complete.  Con- 
sequently, it  is  not  possible  to  expand  an  arbitrary  input  field  distri- 
bution in  terms  of  them,  since  we  can  always  find  field  distributions 
which  are  orthogonal  to  the  space  configurations  of  all  the  discrete 
modes  and  hence  have  vanishing  expansion  coefficients.   It  is  generally 
assiMied  that  this  effect  is  avoided  theoretically  by  making  the  beam 
small  and  neglecting  radial  variation  of  the  fields.  However,  at  least 
one  component  of  the  fields  must  fall  to  zero  at  the  helix,  while  other 
components  are  strong  near  the  helix,  so  the  approximation  may  be  ex- 
pected to  be  considerably  in  error.  Since  there  are  no  other  discrete 
modes,  to  obtain  a  complete  system  we  shall  include  the  continuous 
spectrum  of  the  helical  sheath  system.  This  continuous  spectrum  does 
not  appear  to  have  been  discussed  previously. 

1.  J.  R.  Pierce,  "Theory  of  the  Beam  Traveling  Wave  Tube*,  P.I.R.E. 
_25,  111  (19^'?). 

2.  L.  J.  Chu  ?nd  J.  D.  Jackson,  "IField  Theory  of  Traveling  Wave  Tubes"  , 
P.I.P..E.,  36  ,  853  (19^8). 

3.  J.  R.  Pierce,  "Fffect  of  Passive  Modes  in  TOT 's  ■,  P.I.R.E.,  ^  ,993, 
(19^8). 
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A  second  diffic-alty  arises  from  the  assimption  that  the  fields 
are  circularly  syiimetrical.  The  traveling  wave  tube  is  generally  excited 
by  hPviDg  the  beam  cross  a  waveguide,  and  the  fields  of  the  waveguide  are 
not  circularly  symmetrical  in  the  system  of  cylindrical  coordinates  based 
on  the  beam.  Hence,  components  of  field  having  angulsr  dependence  rnay  be 
generated.  This  possibility  is  included  in  the  analysis  to  be  presented. 

The  method  to  be  used  is  founded  on  the  premise  that  the  inci- 
dent fields  in  the  inpat  waveguide  are  strong  compared  to  the  fields  gener- 
ated by  the  beam  in  the  waveguide.  Consequently,  we  sh^ll  neglect  the  re- 
action of  the  beam  upon  the  input  fields.  As  a  result,  we  may  assume  that 
the  fields  and  the  electron  velocity  and  current  are  given  functions  of 
the  transverse  coordinates  at  the  place  where  the  beam  enters  the  heli- 
cal section.  We  also  assume  that  the  helix  is  sufficiently  long  that  the 
backward  waves  generated  at  the  output  section  may  be  neglected  near  the 
input.  Thus,  we  shall  treat  the  helix  as  semi-infinite  in  length,  with 
the  input  fields  specified  at  z  =  0.  The  question  to  be  treated  is: 
given  the  input  fields,  to  find  the  fields  along  the  helix  and  especi- 
ally to  find  the  strength  of  the  amplified  wave. 

A  natural  method  for  treating  semi-infinite  systems  is  to  apply 
the  Laplace  transformation.  This  approach  will  bring  in  the  initial  con- 
ditions in  a  direct  manner,   ■'■'he  fields  will  be  expressed  as  contour  inte- 
grals, whose  integrands  will  possess  both  poles  and  branch  points.  The 
residues  at  the  poles  will  represent  the  discrete  modes  of  the  system 
and  the  integrals  around  the  necessary  bramch  cuts  will  provide  th?  con- 
tinuous spectrum.  The  theory  may  be  developed  for  arbitrary  angular  de- 
pendence with  no  greater  complications  than  arise  in  the  symmetric  case. 

We  shall  determine  the  current  in  the  beam  at  large  distances 
down  the  helix.  Since  the  radius  of  the  hole  in  the  wall  of  the  output 
waveguide  through  which  the  beam  passes  is  a  small  fraction  of  the  wave- 
length in  the  guide,  most  of  the  field  in  the  helical  section  is  reflected, 
and  the  major  coupling  to  the  output  guide  is  by  means  of  the  a.c.  beam 
current.  This  current  is  expressed  in  terms  of  transfer  coefficients  which 
represent  the  current  that  would  be  generated  by  unit  field,  unit  charge, 
or  unit  a.c.  electron  velocity. 
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In  the  next  section  we  transform  and  solve  the  field  equations, 
The  niatching  conditions  at  the  heam  and  at  the  helix  are  applied  to  the 
solutions  In  Section  3i  where  we  eliminste  the  ar'bitrary  constants  that 
appear.   In  Section  U  we  invert  the  transforms  snd  interpret  the  results. 
In  this  paper  the  analytical  results  are  presented.  Application  to  prac- 
tical exciting  systems  and  some  computations  of  experiiientally  determin- 
ahle  quantities  will  appear  in  a  second  paper. 
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2.  Solvitlon  of  the  Field  Equations 

In  order  to  apply  laplace  transform  theory  we  wish  to  have  the  prop- 
erties of  the  system  independent  of  the  z  coordinate.   Consequently,  we  re- 
place the  physical  helix  hy  a  lossless  cylindrical  sheath  of  radius  a  which 
is  conducting  in  the  helical  direction  only.  This  epproxin>ation  has  "been 
custoEierily  employed  in  traveling  wave  tuhe  research.  We  assume  that  a 
strong  longitudinal  static  magnetic  field  constrains  the  electrons  to  have 
only  longitudinal  motion.  The  lifting  of  this  assumption  will  he  discussed 
in  a  future  y^per.  The  electron  "beam  is  assumed  to  "be  cylindrical,  of  radius 
"b,  and  the  d.c.  electron  velocity  is  ass-umed  to  "be  constant  across  the  "beam. 
This  approxiuvation  limits  the  d.c.  "beam  current,  else  the  d.c.  potential 
will  fall  near  the  "beam  axis.  We  limit  ourselves  to  a  small-signal  analysis, 
and  therefore  assume  that  a.c.  quantities  are  small  compared  to  the  associ- 
ated d.c.  qriantities.   This  approximation  linearizes  the  equations  and  makes 
a  solution  possi'ble.  We  use  cylindrical  coordinates  with  the  z-axis  along 
the  "beam  axis. 

The  system  is  governed  "by  nine  eqtiations;  that  is,  six  Maxwell 
equations,  the  equation  of  continuity,  the  rele.tion  between  current,  charge, 
and  velocity,  and  the  force  equation  for  the  electrons.  We  separate  d.c. and 
a.c.  quantities,  using  the  time  factor  e    for  the  latter,  and  keep  only 
first  order  terms  in  the  a.c.  quantities.  The  equations  "become: 

I       I 

2,1)  7xE,  +  JuiMkg,  =  0 

t        I    I 

~*  I      I 

V  •  £  +  ju»p   =0 

t  t 

m(juj  ■•■  ■''o  37^Z»  =  ^ 
All  primed  quantities  are  a.c,  while  all  quantities  with  subscript 
zero  are  d.c.  We  have  assumed  the  d.c.  electron  velocity  to  "be  much  less 
than  the  velocity  of  light,  and  hence  have  neglected  the  interaction  of  the 
electrons  with  the  a.c.  riagnetic  field.  Because  of  the  strong  d.c.  magnetic 
field,  the  velocity  has  only  a  z  component,  and  the  same  is  true  of  the  cur- 
rent. We  may  save  much  writing  "by  introducing  a  system  of  reduced  units. 
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I  I  I 

We  thus  measiire  v     in  units  of  v   ,    r  in  units  of  r   ,   J     in  units  of 
-J  o     '  ^  0     -»• 

J   =  r  V   .  H     in  units  of  2V   ,  where  V     is  the  heam  voltage,  and  E     in 
o        o   o      — »  0  0  "-* 


units  of  2V    Jtlu..     Thus 


Jtl\i.. 


0 

I  I  I 


V=VV       ,  P=tCP,  J=pvJ 

e/  «  27  E,  ,  H     =  27    v/cTm^  H  . 


»  z  0   z    '  r        ro''    '  'z       ro  o      z    • 


V^e  introduce  some  p.ssocisted  symhols : 

P^^j  k  =  u,/c,  k     =  uj/v^   ,  3  =  v^/c   , 

h^  =  r«/2e7„  =  J  /2€v  7^  =  J„/2e  y2e/m     7^^     . 

'  O     0     0     0  0     0     '         0 

^  .'  -■-  -v.  '- 

The  quantity  h  is  the  dehunching  wave ^familiar  from  klystron  theory.   It  is 
a  measure  of  the  strength  of  the  a.c.  space  charge  forces.  The  assumption 
that  the  d.c.  "beam  velocity  is  constant  across  the  beam  implies  hb  is  small 
compered  to  \mity.   In  these  reduced  units,  the  system  of  equations  becomes; 

V»f^+  jkH^=  0 

7i«H  -  jkE  =  e8i^(p  +  v^)  i 

J  =  p  +  ▼ 

inS 
We  assume  that  all  quantities  have  the  angular  dependence  e 

The  results  may  be  sumr^.ed  over  n  at  the  end  of  the  analysis.  The  next  step 

is  to  take  the  Laplace  transform  with  respect  to  z  of  the  set  of  equations 

2.^).  For  any  quantity  7,   we  define  the  Laplace  transform  T   by: 

2.5)         F(r,e,s)  =  ^F(r,«,2)e"^*  dz 
The  transform  of  a  derivative  is  given  by: 

2.^)  fhY(r,Q,2)   ^-8z  ^^  ^  .F(r,©,0)  +  8F(r.©.8) 
^     oz 

0 


-  6  - 


We  shall  find  the  solutions  of  the  transformed  equations.  We  then  interpret 
the  results  by  using  the  inverse  transform: 


a+ioo 

sz  , 
ds 


2.7)  T(r.e,z)  =  2^  7    ^(^^.^.b)  «' 

c-ioo 
The  contour  lies  to  the  right  of  all  singularities  of  the  integrand. 

The  result  of  transforming  the  equations  Z.^)   is  the  more  manage- 
able set  of  equations: 

^      _ 

2.8)  ^f  -8H^.J1^^  =  .H^(0) 

an 

sH^-  ^  -  JlcE^  =  H^(0) 

-l-rH^-  ^H  -JkS  =  ph^Cp  +  V  ) 
r  ar    Q         r   r   '^  z      '    z 

3v„  +  (s  +  jk')p  =  p(0)  +  v^(0) 

J   =  p  +  V 

(s  +  Jk  )  ?^  =32+  v^(0) 
z      z     z 

The  arguiTient  0  refers  to  z  =  0,  while  all  harred  qxiajitities  are  functions 

InO 
of  r  and  s,  with  the  Q  dependence  e    implied.  We  see  that  the  initial 

conditions  have  come  into  the  eqiaations  in  a  perfectly  natural  manner, 
without  the  necessity  of  matching  them  to  a  modal  decomposition  of  the 
fields. 

The  first  four  equations  of  2.8)  may  be  solved  for  the  trans- 
verse components,  E  ,  E^,  H  ,  H^,  in  terms  of  the  longitudinal  components 

Z^      9      X^      9 

E  ,  H  ,  and  the  incident  field  components.   Similarly,  the  last  three  equa- 
z   z 

tions  may  be  solved  for  the  beam  quantities  v  ,  p,  and  J  in  terms  of  E 

z  z  z 
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and  the  incident  electronic  components.  The  solutions  of  these  seven  equa- 
tions are: 


E  = 


r    2^,  2 
8  +k 


aE 

8  __2.  +  M  H  +  sE  (0)  -  Jk  H.(0) 


3H 


8  +k 


^[^^^J'^  ^*sE^(0)*Jkfl^(0)] 


H 


2.9) 


:  =— i—  -^E  +8  ^+  sH  (0)  +  jkB_(0)  1 
r    2,2  [_r   z      3r     r      *'0*'J 


8  +k 


3E 


^0  =  72^  [-^^  gr  -^  '^  V  «"«^°^  -  J^r^Q)] 


8 


T   =  — 
Z    8 


'^''  (6+Jk') 


try  L- 

+.ik')2  L   2 


P  =  r-^  f-sE^  +  Jk  v^(0)  +  (s  +  jk  )  r(0) 


] 


(8  +2Jk  )  v^(0)  +  (s+Jk  )  p( 


.,] 


The  solutiofts  2.9)  may  "be  inserted  into  the  two  remelning  equations  of  2.8). 
We  define  the  "radial  propagation  constants  "  p  and  \    "by  the  relations: 


2.10) 


2     /  2  ^  ,  2v 
p  =  _  (s  +  k  ) 


2    2 

\  =  P 


1  + 


(s 


'j.')!) 


The  terms  involving  the  incident  fields  may  he  considerably  simplified  hy 
tisin,?-  the  Maxwell's  equation  2.^).  The  transverse  input  fields  may  be 
elimin8.ted,  and  the  equations  expressed  in  terms  of  the  longitudinal  in- 
put fields  snd  their  derivatives.  The  resulting  wave  equations  ere: 

2- 


r 


'dz 


*-  r  -•  Jk  *-    (a  +  Jk  ) 

r(0)]  . 


'2,2 

k  -k 

8  +  Jk' 
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The  second  of  these  equations  holds  only  for  r  <  1).  V-Tien  t  >  Ta ,  so  v;e  are 
outside  the  heam,  the  input  velocity  end  charge  density  are  zero,  and  ->,_  is 
replaced  hy  p. 

We  now  wish  to  solve  the  wave  equations  2.11)  regarding  the  right- 
hand-sides  as  forcing  terms,  which  are  given  functions  of  position.  Bound- 
ary conditions  are  applied  at  the  "beam  radius  r  =  b  and  the  sheath  radius 
r  =  a.  The  exact  "boundary  conditions  to  he  used  will  "be  stated  in  Section 
3.  To  solve  the  equations,  we  introduce  the  Green's  functions: 


r 

< 

1 
r 

r 

> 

1 
r 

r 

< 

r 

r 

> 

r 

^1=  V^^^^n^^^  ) 

K  (pr)  I  (pr') 
n     n 

We  are  assuaing  that  K  (pr)  is  that  solution  of  the  homogeneous  wave  equa- 
tion which  is  exponentially  decreasing  at  infinity.  This  requires  that  p 
have  a  positive  real  part.  Howeve.-,  p  is  defined  in  terms  of  the  coarplex 
variable  s,  and  s  lies  on  a  vertical  contour  in  the  right  half -plane.  As 
a  function  of  s,  p  hs.s  branch  points  at  s  =  jk,  s  =  -jk.  The  branch  cuts 
may  be  so  placed  that  p  actually  does  have  a  positive  or  zero  real  part 
along  the  entiire  path  of  integration  in  the  s-plane.  The  details  of  this 
analysis  will  be  presented  in  Section  k. 

We  denote  the  right-hand  sides  of  Equations  2.11)  by  7^  snd  I'j.+I'j. 
respectively.  Applying  Green's  theorem  to  the  f\inctions  E^,H^,  and  the 
appropriate  Green's  functions  yields  the  solutions  of  2.11),  each  in  its 
appropriate  region  of  space,     -^  — 

^zi  '^   ^^n^'^^^  ~^  r'dr'G^(r,r')[l'j,(r')  +  TjCr')]      r<b 

^zm  '^  ^2^n^^^^  *   AoK^Cpr)-/  r'dr'G^Cr.r')  ?j,(r')     b  <  r  <  a 

b 
2  13)  00 

i   =  A^K^(pr)-/  r'dr'GgCr.r')  TjjCr')  r>a 


^.i.m  =  B,I^(pr)-/r'dr'G2(r.r.)r„(r') 

0 

H   =  3,K  (pr)-/°°r'dr'G-(r.r')  F  (r') 
zo    <i  n     'L       c  n 


r  <  a 


r  >  a 
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The  letters  i,m  nnd  o,  stand  for  inner  region  r  <  "b,  middle  region  "b  <  r  <  a, 

and  outer  region  r  >  a,  respectively.  Since  the  longitudinal  magnetic  field 

is  unaffected  "by  the  presence  of  the  "beam,  we  may  tise  the  same  solution  for 

H  in  regions  i  and  m.  The  integrals  appearing  here  represent  particular 
z 

integrals  of  the  equations  2,11)  due  to  the  forcing  terms,  while  we  add 
appropriate  solutions  of  the  homogeneous  equations  to  satisfy  the  "boundary 
conditions.  The  conditions  of  regij.larity  at  r  =  0  and  r  =  co  are  already 
satisfied  hy  the  solutions  2.13) •  We  shall  now  use  the  "boundary  conditions 
at  the  "beam  radius  and  the  sheath  radius  to  eliminate  the  arbitrary  con- 
stants A^,...B-.  This  reduction  will  form  the  su"bJeot  matter  of  the  next 
section. 
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3.  Application  of  the  Boundary  Conditions 

The  ■botmdary  conditions  at  the  beam  radius  are  that  the  tangen- 
tial electric  and  magnetic  fields  shall  he  continuous.  The  second  condi- 
tion Is  a  conseq.uence  of  the  presence  of  the  infinitely  strong  static  mag- 
netic field.   If  the  static  field  were  finite,  radial  spreading  of  the 
"beam  woiild  occur,  and  this  would  have  to  he  replaced  hy  a  surface  current 
density  at  the  "beam   surface.  This  surface  current  would  cause  a  discon- 
tinuity in  the  tangential  magnetic  field.  However,  this  effect  is  sup- 
pressed "by  the  infinite  static  field.  Since  the  boundary  conditions  are 
to  he  satisfied  for  all  z  and  all  ©,  the  transfojnned  field  components  mtist 

also  satisfy  the  conditions.  The  continuity  of  H  is  already  insuired  by 

z 

the  solutions  2.13).  Hence,  the  boundary  conditions  to  be  applied  at  the 
beam  are: 

3.1)  zi  ~  zm 

Eo*  =  E^      r  =  b 

^01  "  ^Om 
In  view  of  the  relations  2.9).  these  conditions  can  be  replaced  by  the 
equivalent  set  of  conditions: 

3.2)  \i°\m 

8r   zi    Sr  zm 
At  the  helical  sheath,  the  boundary  conditions  are  that  the  com- 
ponent of  electric  field  in  the  helical  direction  shall  vanish,  the  component 
of  electric  field  upon  the  sheath  but  normal  to  the  helical  direction  shall 
be  continuous,  and  the  component  of  magnetic  field  in  the  helical  direction 
shall  be  continuous.  The  transformed  components  nrast  also  satisfy  these 
conditions.  The  helix  direction  makes  the  angle  tt/2  -  a  with  the  z-axis, 
whence  these  conditions  are  expressed  as: 

E-.  cos  a  +  B   sin  a  =  0 
jOm        _zm 

-  ^v     Ert  cos  a  +  E   sin  a  =  0 

3.3)  Oo         zo  ^  ^  ^ 

E-_  sin  a  -  E   cos  a  =  E.  sin  a  -  E   cos  a 
9m  zm        Qo  zo 

E-.     cos  a  *  H   sin  a  =  H„  cos  a  *■  H   sin  a 
dm        zm        do         zo 
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Upon  using  the  relations  2.9) i  these  conditions  may  be  simplified  slightly, 

yielding: 

i   =  i 
zm    zo 

3.4)  Sr   zffl  ~  Sr  zo 

B-     cos  a+B       sina=0 
Om  zm 

H-  cos  a  +  H   sin  a  =  H„  cos  a  +  H   sin  a 
dm         zm         Qo  zo 

Thus  equations  3.2)  and  J.k)   provide  six  conditions  to  "be  satisfied  by  the 
solutions  2.13)-  These  solutions  contain  six  arbitrary  constants,  and 
since  the  forcing  terms  render  the  boundary  condition  eqtiations  inhomo- 
geneous,  a  solution  can  be  found  for  the  arbitrary  constants  if  the  deter- 
minant of  the  equations  does  not  vanish.  Since  the  determinant  will  in- 
volve Bessel  functions  of  the  complex  variable  p,  we  may  expect  that  the 
determinant  will  vanish  only  for  certain  values  of  p,  if  any.  These 
critical  values  will  correspond  to  the  discrete  modes  of  the  system,  which 
can  be  set  up  without  forcing  terms.  Thus,  the  determinant  does  not  van- 
ish identically. 

The  boundary  condition  equations  may  be  simplified  by  introduc- 
ing a  set  of  "transfer  coefficients".  These  refer  to  the  values  of  the 
integrals  appearing  in  2.13)  when  r  =  a  or  b,  and  will  represent  the 
effect  of  the  initial  conditions.  The  six  coefficients  we  introduce  are 
defined  by: 


^l=/^<^^^n(n^^  [^E^^^*^J<^0 


0 


3.5) 


C^  =y  r  dr  K^(pr)  Fj.(r) 

b 
S=/^<^^  I>r)Pj.(r) 

D  .00 


C^=  /  T  AT  Kjvr)  l^ir) 

a 

^1=/'^^^  I^(pr)r„(r) 

0 

D2=/rdrZ^(pr)rj^(r) 
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The  integrrals  that  appear  here  will  all  converge  if  J-  and  r„  are  suffi- 
ciently regular  at  infinity,  and  this  condition  is  always  satisfied  for 
practical  exciting  systems.  We  also  introduce  a  constant  |ji,  defined  ty: 

3.6)  M-  =  P  tan  a  -  jn  a/a  . 

The  six  "boundary  condition  equations  niay  now  "be  written  as: 

'iM^I^(^b)-pM2l^(pl))-pM^^(pU)  =\-b  C^E^(^T))-  pb  C^  I^(pb) 

Aglj^Cpa)  +  A^Kj^(pa)  -  A^K^(pa)  =  C^^(pa)  -  C^I^(pa) 
B^I^(pa)  -  B^I^^Cpa)   =  D^^K^Cpa)  -  ^2^'j^^^ 

3.7) 

-tilA^Ij^Cpa)*  A^^(pa)~|  +jkpB^I^(pa)  =  -^iC^^(pa)+JkpD^K^(pa)-(sE^(0,a)+jkHj.(0,a)} 
jkpU^I^Cpa)  +  A^^(pa)-A^K^(pa)  J  +  ^ik  I^(pa)-B2K^(pa)J     = 

jkp|c^^(pa)-  C^I^(pa)J  +  ^   [^i^n^^^  ~  ^2^n^^J 

The  prime  denotes  the  derivative  of  a  Bessel  function  with  respect  to  its 
argument . 

The  qiiantity  which  is  to  be  determined  is  the  current  in  the  beam. 
From  the  solutions  2.9) i  the  transformed  current  in  the  beam  is  propor- 
tional to  the  transformed  longitudinal  electric  field  E  .  Consequently, 

z 

the  constant  which  is  to  be  solved  for  is  A^,  which  measures  the  strength 
of  the  electric  field  in  the  beam.  However,  it  proves  most  convenient  to 
solve  equations  3.7)  for  A-,  and  then  express  A,  in  terms  of  A^.  We  shall 
proceed  to  solve  equations  3.7)  in  sequence,  rather  than  by  using  determin- 
ants, since  for  a  sixth-order  system  such  as  3.7)  the  determinants  become 
too  unwieldy. 

First,  we  solve  the  third  equation  of  3.7)  for  A^,  and  then 
solve  the  fourth  eqiaation  for  B-.  This  expresses  the  constants  describ- 
ing the  fields  outside  the  helical  sheath  in  terms  of  the  constants  for 

for  the  region  inside  the  sheath.  There  results: 

I  (pa) 

3.8)  A^  =  (A^*  C^)  YJ^    *  (^3  -  °3^ 

n  " 

\   (pa) 
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Next,  these  two  relations  are  inserted  into  the  last  of  equptions  3.7). 
The  resulting  equation  nay  be  solved  for  B^  in  terms  of  A_  and  A„.  We 
thereby  obtain  the  magnetic  field  constants  in  terms  of  the  electric 
field  constants.  The  result  may  be  reduced  by  the  Vfronskian  relation- 
ship: 

^•^^  i'(0  kU)  -  lU)   k'(z)  =  1/z 

n     n      n     n 

and  the  solution  becomes: 

41    K  (pa) 

This  may  be  reinserted  into  3»^)  to  express  the  outer  niagnetic  constant 
B-  in  terms  of  electric  field  ccnstento,  yielding: 

,,    I '(pa) 

Thus,  to  this  point  the  outer  electric  field  and  the  entire  magnetic  field 
ha.ve  been  eliminated,  and  we  are  left  with  three  equations  for  the  inner 
and  middle  electric  field. 

The  third  step  in  the  reduction  is  to  eliminate  the  constants 
A.  and  A_  by  solving  the  first  and  the  second  of  equations  3«7).  We 
introduce  the  symbol: 

3.1?)        6  =\b  I^(Yb)  K^(pb)  -  pb  Ij^(\T5)  K^(pb) 

This  quantity  6  is  the  deteminant  of  the  first  two  equations  of  3.7).   It 
reduces  to  unity  when  the  beam  current  goes  to  zero,  so  that  h  =  0  andii=  p, 
The  solutions  for  A,  and  A.  are  found  to  be': 

3.13)      h  =  l  W'^Z  "^  '^l^'^^  K^('\^)Kj^rpb)  -  pb  K^(pb)  K^('.ib)>  J 

A3  =  I"  [-Cj^  +  (A2-C2){pb  I^(pb)l^(->ib)--Yb  I^(\^)lj,(pT3)}] 

It  is  essiimed  here  that  6  does  not  vanish.   Since  for  large  values  of  Res 
the  quantity  "Y  as  defined  by  2.10)  is  equal  to  p  +  0(l/s),  for  such  values 
6=1+  0(1/8).   Consequently  we  may  select  a  contour  of  integration  such 
that  6  does  not  vanish  anywhere  along  it. 
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The  penultimate  step  is  to  take  the  relations  3.8) i  3.10),  3.11) • 
pnd  3.13);  and  substitute  all  of  them  into  the  fifth  of  equations  3.7). 
The  only  remaining  undetermined  constant  is  A  .  This  resulting  equation 
may  be  simplified  at  nfi-ny  points  by  use  of  the  Wronskian  relation  3.9) » 
and  several  other  terms  cancel.  We  define  the  function: 

I^(pa)K^(pa)  +  i^£  I^(pa)K^(pa) 

3.li+) 


Gj^(pa)  = 


K2(pa) 


The  solution  for  A  may  be  expressed  in  the  form; 


3.15) 


Ag  =  N2/M2 


"2  "  ^n^^^  +{pT5lj^(pb)  lj\-b)   -^bl^(^b)  Ij^(pb)}  /6 
^2"  ^1/^  *  C2{pbl^(pb)  I^(iib)-^bl^(-'ib)  I^(pb))  /6 
+  C^+  (ljj(pa)-G^(pa)K^(pa)]  C^,\(pa) 

K'(-pa)        1    I  (p^)         1 


K  (pa) 


K 


The  final  step  is  to  substitute  this  result  into  the  first  of 
eqijations  3.13)  to  obtain  A^.  Only  the  coefficients  of  C.  and  C^  In  N_ 

will  be  affected.  After  some  reduction,  the  final  result  for  A^  is: 

3.16)  A^  =  N^/M^ 

M^  =  pb  Ij^(^b)  {  I^(pb)  -  G^(pa)  K^Cpb)} 

-  lib  I^(\b){  I^(pb)  -  G^(pa)  K^(pb)) 

pb  \{\^)   (I^CP^)  -  G^Cpa)  Kn(P^)}] 

°1 
-\bK^(^b){l^(pb)  -  G^(pa)  K^(pb)) 

-  G^(pa)C2  +  C^  +  {In(pa)  -  G^^Cpa)  \^V^))   <^J\^y^) 
,   K  (pa)  I  (pa)       T 

n  n  n 


^1  = 
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With  this  determine tion  of  the  constant  A^,  the  transformed  electric  field 

E  is  known  completely.  The  last  of  equations  2.9)  now  says  that  the  trans- 

^  —II? 

formed  ciirrent  in  the  beam  is  found  by  multiplying  E  by  Jk  /(8+  Jk  )  and 

adding  some  initial  terms.  The  next  step  is  to  invert  the  Laplace  trans- 
formation to  obtain  the  actual  current  and  then  interpret  the  results. 
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I|.  Besulta  and  Interpretations 

We  now  shall  invert  the  Laplace  transformation.  The  first 
steps  in  this  process  are  to  properly  locate  in  the  complex  s-plane  the 
contour  of  integration  and  any  "branch  cuts  that  may  he  required.  The  con- 
tour DTUSt  he  placed  to  the  right  of  all  sin^nilarities,  so  we  mist  study 
where  the  singularities  are  located. 

The  first  step  is  to  demonstrate  that  we  can  actually  satisfy 
the  requirement  that  the  real  part  of  p  is  to  he  positive.  The  defini- 
tion p  =  -(b  +k  )  shows  that  p  has  branch  points  at  s  =  i  Jk.  Since 
the  contour  of  integration  is  to  lie  in  the  right  half-plane,  the  hranch 
cuts  which  rrnist  he  inserted  to  render  p  single  valued  cannot  extend  to 
+  00,  We  write  s  =crt  it  .  By  its  definition,  p  is  purely  imaginary  on 
the  entire  real  axis.t  =0.   It  is  also  purely  imaginary  on  the  imaginary 
axis(r=0  when  ]  r  |<  k,  and  is  purely  real  on  the  imaginary  axis  when  \r   |>  k. 
The  requirement  is  that  |arg  p|  ^  tt/2,  so  that  the  K  function  either  is 
exponentially  small  or  an  outgoing  wave  when  r  approaches  infinity: 

We  shall  place  the  branch  cuts  as  in  Fii'ure  1. 


1 


T 


s-  plane 


-jk 


The  Branch  Cuts  for  p 
Figure  1 

The  upper  hranch  cut  starts  from  s  =  Jk,  goes  along  the  imaginary  axis 
to  zero,  remaining  in  the  second  quadrant,  and  then  passes  to  minus  in- 
finity along  the  upper  side  of  the  negative  real  axis.  The  lower  hranch 
cut  Is  a  mirror  Imare  of  the  upper,  and  it  lies  wholly  in  the  third  quad- 
rant. 


-  17  - 


We  now  define  the  square  root  involved  in  p  "by: 

s    TT       

2    /  2     2 
p  =  e        \/s+k  't  =  lins>0 


^  y~2    2~ 

=  e  /s  +  k        t=lms<0 

/"i   2 
The  square  root  appearing  here  is  so  taken  that  lim  i/s  +  k  =  k,  and  it 

is  real  positive  along  hoth  branch  cuts.  Hith  this  definition,  p  is  posi- 
tive imspinary  along  the  upper  branch  cut,  has  a  phase  between  zero  and  v/Z 
in  the  second  quadrant,  and  is  positive  real  along  the  upper  section  r  >  k 
of  the  positive  imaginary  axis,  '^he  phase  of  p  lies  between  zero  and  -tt/2 
in  the  first  quadrant,  and  p  is  negative  imaginary  on  the  upper  side  of  the 
positive  real  axis. 

In  the  lower  lialf  plane,  p  is  negative  imagine.ry  on  the  lower  branch 
cut.   Its  phase  is  betv;een  zero  and  -tt/2,  in  the  third  quadrant,  and  p  is  posi- 
tive real  along  the  lower  section  -t>  k  of  the  negative  iine.gina.ry  axis. 
In  the  fourth  quadrant,  the  phase  of  p  is  between  zero  and  tt/2,  rnd  p  is  posi- 
tive imaginary  on  the  lower  side  of  the  positive  real  axis.  There  accordingly 
is  p   jtunp  in  the  value  of  p  upon  crossing  the  real  axis.  However,  this  Jun^) 
does  not  contribute  to  the  value  of  an  integral.  We  see  that  with  the  branch 
cuts  of  Figure  1  and  the  definition  ^.l)  we  actually  have  satisfied  the  re- 
quirement jarg  p|  S  ""/Z  . 

Other  branch  points  of  the  integrand  occur  at  the  zeroes  of \   . 
By  the  definition  2.10),  \  has  branch  points  at  3  =  +  Jk  arising  from  the 
zeroes  of  p.   It  also  has  branch  points  at  s  =  -j(k  ^  h),  and  a  simple  pole 
at  s  =  -jk  .  V/e  shall  dispose  of  these  latter  singulrrities  by  a  closed  con- 
tour which  surrounds  the  three  singular  points.  The  phase  of  the  square 

2 

root  of  the  second  factor  in  \     is  specified  by  requiring  that 

r     2/        '  2"i^/^ 

1  +  h  /(s  +  Jk  )      approaches  1  uniformly  with  respect  to  arg  s  as  Isj 
approaches  infinity.  This  condition  may  be  satisfied  by  the  closed  contour 
branch  cut.  Accordingly,  the  path  of  integration,  L,  and  the  branch  cuts 
are  given  as  in  Figure  2 
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4 

kT 

i 

s-plone 

1 
< ^ ' 

jk 

0 

"«             r 

<r 

V 

-jk 

^-j(k'-h) 

•  -jk 
J-j(k'+h) 

\                     u 

The  Branch  Cuts  for  \ 
Figure  2 

The  integrand  has  hranch  points  at  -J  (k  -h)  -nd  -j(k+h),  end  apparently 

t 
an  essentisl  singularity  at  -Jk  .   Consequently,  the  indicated  section 

of  the  imaginary  axis  must  be  excluded  from  the  complex  s-plene. 

We  would  like  to  evaluete  the  Laplace  inversion  integral  2,?) 
by  closing  the  contour  in  the  left  half  plane  and  finding  the  contribu- 
tion from  the  various  singularities.  The  contour  may  only  be  closed  If 
the  integral  converges  over  a  large  semicircle  in  the  left  half  plane. 
Accordingly,  the  next  step  is  to  investigate  the  behavior  of  the  inte- 
grand on  this  semicircle. 
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We  shall  first  dispose  of  the  second  term  in  eqijation  2.13)  for 

B  .  This  contributes  to  J  the  amotint: 
z  z 

..2)  J,,  =  -/r'dr'[r,(r')  *  .,(r')].  ^/ -Ji^  ,"  I„(^r')V,r) 

-/Vdr'[r,(r')  *  .,(r')]  ^  /  f^^  '"  ^(^')  ^.(V'") 

The  two  integrals  over  s  are  of  the  same  character,  since  in  each  the  argu- 
ment of  the  I  function  is  smaller  than  the  argument  of  the  Z  function, 
n  1 

On  the  large  semicircle  we  have'\=  p  +  oCl/s),  and  Re  p  >  0.  Consequently, 
on  the  upper  quartercircle  we  apply  the  asymptotic  forms: 

h.3)  e   I^(^r  )K^(^r)  ->e     ^=;:  '  7^  e  ^ 

1      8z-p'(r-r') 


2p/rr' 

Since  Re  p  >  0,  Re  p(r-r')>0,  for  r>r'.  Also,  Re  s  <  0.  Therefore,  the 
integrand  is  exponentially  small  on  the  upper  quartercircle.  The  same 
reasoning  is  valid  for  the  lower  quartercircle.  Also,  the  integrals 
along  straight  lines  connecting  the  imaginary  axis  with  the  contour  of 
integration  will  also  "be  exponentially  small  as  the  lines  move  to  infin- 
ity, except  Rt  the  imaginsry  axis  itself.  At  the  axis,  the  integrand  is  of 
order  s"^^,  and  there  is  no  contrihution.  Accordingly,  the  integral  for  J 

may  he  closed  in  the  left  half -plane. 

The  treatment  of  the  A,  I  (^r)  term  is  more  complicated.  For 

large  p,  the  denominator  M,  in  A^  reduces  to  G  (pa).  The  asymptotic  form 

for  G  (pa)  is: 
n  -^ 

U,U)  G^(p8)  ->  ie^* 

The  several  terns  in  the  numerator  N,  shall  be  treated  separately,  and  we 

-2pa 
shall  show  that  none  of  them  is  large  enough  to  cancel  the  factor  e  '^ 

which  arises  from  1/M^. 
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We  first  discusB  the  order  of  magnitude  of  the  C,D  coefficients. 
From  the  definitions  2.5)  we  obtain,  vcpon  using  the  monotonic  properties 
of  the  I  and  K  functions  for  |p|  sufficiently  large  and  Re  p  >  0: 

k.5)   C^  -^e^^  .  C^  ^e-P\  C3->e^.  C^-^e'^^   .  D^-^  e^*  .  D^  -^e'^^ 

The  bracketed  factor  multiplying  C^^  is  of  the  order  1  +  ot^G  (pa)e"'^  ,  where 

CL.    is  an  algebraic  function  of  p.  Multiplying  this  by  C,  gives  eaqjressions 
of  the  order  e   and  e*^   .  Both  of  these  lead  to  exponentially  small 
terms  when  divided  by  M, .  Denoting  the  remaining  terms  in  A,  by  R^,  R-, 
R^,,  S.  ,3-,  we  may  show  directly  that  R^  is  of  the  order  e~^  ,  and  the 

other  terms  are  of  the  order  e  ^^ .     Therefore,  all  the  terms  in  A^  are  ex- 
ponentially small  on  the  semicircle  and  the  integral  may  be  closed  in  the 
left  h^lf-plane  and  replaced  by  an  integral  around  the  singularities  of  the 
integrand. 

The  only  possible  singularities  of  the  integral  are  branch  points, 
which  are  taken  care  of  by  the  branch  cuts  of  Figure  2,  a  possible  essential 
singularity  at  s  =  -Jk  ,  which  we  shall  discuss  later,  and  poles  at  the  zeroes, 
if  any,  of  M, .  The  zeroes  of  M,  correspond  to  the  discrete  modes  of  the  sys- 

1  1        ^  B^z 

tern,  since  evaluating  the  residues  leads  to  terms  behaving  like  e   ,  where 
s.  is  the  value  of  s  at  the  i'th  pole  p.  of  the  integrand. 

This  shows  the  great  advantage  of  the  Laplace  transform  method. 
The  amplitudes  of  the  discrete  modes  are  given  by  the  residues  of  the  inte- 
gral at  its  poles.  No  matching  of  solutions  in  the  helical  system  to  input 
fields  need  be  employed.   Instead,  the  coinplete  description  of  the  matching 
is  given  by  the  coefficients  C^,...D  ,  which  are  effectively  the  Fourier 
coefficients  of  the  input  fields  with  respect  to  the  modes  of  the  helical 
system.  Since  these  modes  are  not  orthogonal,  the  eu;tual  solution  of  the 

matchinf  equations  would  be  a  very  difficult  task.  A  partial  treatment  has 

Li    5 
been  given  in  papers  by  Friedman  ejid  Malin  '  . 

k,   B.  Friedman  and  H.  Malin, '(Excitation  of  the  Traveling  Wave  Tube",  N.Y.U. 

Mathematics  Research  Group  Report  No.  TW-12,  February,  1950, 
5.  B.  Friedman  and  H.  Malin,  "Orthogonality  Properties  of  the  Modes  of  the 

Helical  Wire  Guide",  N.Y.U.  Mathematics  Research  Group  Report  No.TW-13, 

October,  1950. 
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The  trenscendental  equation  descriting  the  discrete  modes  of  the 
system  is  given  ty  setting  M^  =  0.   It  may  be  written  in  the  form: 

0  =  ph  I^(\b)  n  (pb)  -  G^(pa)  K^(pt)^ 

^•^^  ^^   I^(\b)     I^Cpb)  -  G^(pa)  Z^(pb) 


^(\T3)     IjJCpb)  -  G^(pa)  Z^ 
p  I  (\b)     I  (pb)  -  g' (pa)  K  (pb) 


The  vanishing  of  the  s econd  bracketed  expression  in  ^.6)  leads  to  the  trans- 

2 
cendental  equrtlon  derived  by  Chu  and  Jackson  and  given  in  simplified  form 

as  equation  30)  of  their  paper  in  the  special  case  n  =  0.  The  apparent  modes 

associated  with  the  vanishing  of  I  ("Hb)  or  I  (pb)  -  G  (pa)  K  (pb)  are  spuri- 

n  *-       n       n      n 

ous,  since  for  them  the  second  bracket  becomes  infinite  anc;  M,  does  not 
actually  vanish. 

The  determination  of  the  solutions  of  equation  ^.6),  viewed  as  an 
equation  in  p,  is  a  highly  complex  task.  We  shall  accordingly  use  a  physical 
argument,  which  cannot  be  completely  Justified.   In  the  absence  cf  the  beam, 
equation  ^.6)  reduces  to  G  (pa)  =  0.  The  roots  of  this  simpler  eq^aation  have 
been  discussed  in  greet  detail  by  Phillips  and  Malin  .  We  shall  assume  that 
in  the  presence  of  the  beam,  each  real  root  of  G  (pa)=  0  which  corresponds  to 
forward  propagation  is  split  into  three  roots,  one  real  and  two  complex,  ex- 
actly as  in  the  case  n  =  0.  Since  a  cubic  equation  for  p  may  be  derived  from 
h,6)   by  an  expansion  to  the  first  power  in  pov/ers  of  the  beam  current,  this 
supposition  is  reasonably  well  justified,  although  a  proof  that  the  cubic 
actually  has  complex  roots  has  not  been  given. 

In  the  case  n  =  0,  the  equation  G  (pa)  =  0  has  one  real  positive 
root.  It  has  no  other  roots  in  the  complex  p-plane  for  Ee  p  >  0,  thoiogh  there 
are  roots  for  Ee  p  <  0  in  the  neighborhood  of  pa  =  1  jx  ,  where  x     is  the 
n'th  root  of  J  (x)  =  0.  These  roots  do  not  lie  in  the  allowed  portion  of 
the  p-plane.  Hence,  for  n  =  0,  the  equation  M.  =  0  h?,s  just  three  roots, 
of  which  one  is  real  and  positive  and  the  other  two  are  complex  conjugate. 
Denote  these  by  p^  =  x,  +  iy-if  Pj,  =  x,-  iy^,p^  ~  ^"^  * 

We  now  refer  to  the  location  of  the  poles  in  the  complex  s-plane. 

Corresponding  to  p,  there  will  be  a  pole,  s^  in  the  fourth  quadrant.  The 

value,  s_  of  s  corresponding  to  Pp  is  in  the  third  quadrant.  There  are  also 

_ 

6.  R.  S.  Phillips  and  H.  ?-telin,  "  A  Helical  Wave  Guide  II,  N.Y.U. ,  Mathe- 
matics Research  Group  Report  No.  170-3,  August,  19^7. 
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two  poles  on  the  imaginary  axis  at  the  conjugate  complex  values  s„  ,  8„  . 
The  root  at  s^  gives  an  amplified  forward  wave;   that  at  s  an  attenuated 
forward  wave;   that  at  s-  a  wave  that  propagates  forward  without  amplifi- 
cation or  attenuation;  and  that  at  s-^  a  "backward  propagated  wave.  The 
backward  wave  is  usually  ignored,  but  it  is  required  here  to  satisfy  the 
initial  conditions.   It  might  be  removed  by  providing  lumped  or  distributed 
attenuation,  but  that  is  another  problem. 

'.^en  the  tube  Is  tuned  to  resonance,  the  value  of  a  corresponding 
to  the  root  x  of  G  (x)  =  0  coincides  with  the  singularity  at  8  =  -jk  .  For 

0      0 

small  values  of  h,  the  roots  3^,s„,8p„  of  the  system  with  the  "beam  will  lie 
close  to  this  point  in  the  complex  plane.  The  roots  s^  and  s  ,  because  of 
their  non-vanishing  real  parts,  may  be  placed  outside  the  branch  cut  that 
encircles  the  section  of  the  imaginary  axis  from  s  =  - j (k  -h)  to  s  =  -j(k'+h). 
However,  3_  lies  on  the  imaginary  axis,  and  its  position  with  respect  to  the 
singularity  should  be  investigated.  However,  t}ie  wave  corresponding  to  s  _  is 
a  non-amplified  wave,  and  it  may  be  associated  with  the  bundle  of  non-emplif ied 
waves  that  will  arise  from  the  integration  arotind  the  branch  cut. 

The  amplified  wave  arises  from  the  pole  at  s^.  When  we  evaluate 
the  residue  at  this  pole,  and  simplify  slightly,  there  results: 


i|.7)      J  -      "^ 


8,  N^(8^) 


-  (^Jv^b)   Cg  +  C^  +(k  cot  a)^I^(p^a)K^(p^a)  C^/  p^  K^ (P^a) 
+  jkcotaK^(p^a)  D^/pj^K^(p^a)-Jkcot  a  Ij^(p^a)  C^/Pj^K^CPj^a) 
+  (8Eg(0,a)  +  jkH^(0,a))/  p^  tan  a  K^(p^a) 
M^(8^)  =  SM^(8)/a8  I 

8=8^ 
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Since  s^  is  the  propagation  constant  of  the  amplified  v;ave,  and  M,  is  prop- 
ortional to  the  discontinuity  in  admittance  of  the  system  upon  crossing  the 
"beam  radius,  we  see  that  the  cnarrent  in  the  beam  is  inversely  proportional 
to  the  derivative  of  the  admittance  discontinuity  with  respect  to  propag-a- 
tion  constant,  evaluated  where  this  discontinuity  vanishes. 

The  numerator  N, (s- )  shows  how  the  initial  quantities  induce  cur- 
rent in  the  heam.  We  may  split  N.  into  the  sum  of  three  expressions,  which 
are  respectively  the  contributions  of  initial  currents,  initial  longitudinal 
electric  fields,  and  initial  longitudinal  magnetic  fields.  These  expressions 
are: 

^.8)      ^J  =  {^o^Pl^^-^o^Pl'^^o^Pl^^)/  ^  ^^  ^j(r)^o^^l^^/^o^^l^^ 


'hi   =  {^o^^l^K^Pl^^^oV^^/^  ^^  ^E^^)  ^o^U^^/^o^^l^^ 

0 

a 

+  /r  dr  yj,(r)  k(Pir)-G^(Pia)K^(p^r)] 
"b  *-  -* 

P  T  2   2 

+  (k  cot  a)^  I^(p^a)Z^(p^a)/  r  dr  fj,(r)  Z^(p3^r)/p^  ^o^^l*^ 

a 
-a 

/r  dr  Fj^(r)  I^(p^r)K^(p^a) 
_o 

/  r  <ir  \(^^K^^y''^   IiCPi^)]  /PiKjp^a) 


U,,  =  Jkcot  a 

M 


+ 

a 


+  (3^E^(o,a)  +JkH^(0,a))/pJ  taji  a  K^CPj^a) 


The  expression  for  1I„  may  he  simplified  by  a  judicious  integration 
"by  parts,  which  will  eliminate  the  extra  non-integral  term.  This  partial 
integration  may  he  perform.ed  in  the  expression  for  U,  for  general  n.  The 
resulting  expression  then  attains  a  very  neat  and  cogent  form  if  the  inci- 
dent fields  are  expressed  in  rectangular  components  and  if  the  Q   dependence 
of  the   incident  fields  is  restored.  There  results: 


-  Zh  - 
0   0     L 

+  ^jTofv   dr  t  ;L(pr)e''^^-^^®{s(E^-jEy)+k(H^-jHy)} 
0   a.     ^ 

in  which  the  argpj^ient  of  E  ,  etc.  is  r,9.  Thus,  the  incident  fields  are  split 
into  circularly  polarized  waves.  The  contribution  of  these  circularly  polar- 
ized waves  is  then  found  "by  multiplying  "by  the  appropriate  Green's  function, 
as  specified  iDy  ^.9).  "Proa  the  setup  of  ^.9),  we  see  that  the  amplitude  of 
the  n'th  Fourier  coefficient  of  the  current  in  the  heam  depends  upon  the 
(n+l)'3t  and  (n-l) 'st  components  of  incident  transverse  electric  field. 

A  similar  transforraation  may  "be  performed  upon  IL.   One  of  the  terms 
which  arises  from  the  integration  "by  parts  is  proportional  to  M^^,  and  hence  it 
does  not  contribute  to  the  discrete  spectrum.   It  may  "be  associated  with  the 
integral  4.2)  for  J  „.  The  other  terms  cancel.  Again  restoring  the  Q   depend- 
ence, the  formula  for  Kg,  analogously  to  4.9)  for  Uj^  becomes: 

4.10) 

jd«/r  dr[l^__^(v)e'^^''-^^^{s(E^-JH  )-k(E^-jEy)} 
0   0     '- 

-^%'»fr  ^r[{l„.i(prH(!^(pa)K^.,(pr)]  f"  <''-^"'{sCH^-JHy)-k(VJi:y)} 

0    b    *- 

-{^n*l(P^^*-n^P-^^n*l^^-^'^^"*'^'{^^V^^^-^^'^(V-^y^)] 

-  ^(^n^P^^-^n^P^^^n^P^^-^^^/^  drrz^_^(pr)eJ^^-^^^(s(H^-JHy)-k(E^-jSy)} 

0  a     ^ 

-K^^^(pr)eJ(^*^>«{3(H^+jHj^)+k(E^H-JEy)}-]  /K^(pa). 


-  25  - 

For  the  discrete  spectrum,  the  coefficient  of  the  first  term  simplifies  to 

-"^,(1  (p-b)-G  (p,a)K  (p-jh)}  /gnk  I  (-^^h).  Since  for  the  discrete  spectrum 

8  is  not  far  from  -Jk  and  accordingly  js|  »k  for  slow  waves,  we  are  Justi- 
fied in  referring-  to  the  contrihution  of  N,,  as  predominantly  transverse  electric, 
since  in  the  integrand  s  multiplies  the  transverse  electric  field,  while  in  N 

s  multiplies  the  transverse  magnetic  field  and   !!_  is  predominantly  transverse 
magnetic.  Equation  4.10)  gives  IL,  ns  a  function  of  the  trrnsverse  input  fields. 
The  equivalent  representation  in  terms  of  the  longitudinal  electric  input  field 
is  given  "by  k.8)   in  the  case  n  =  0.  However,  4.10)  is  useful  in  shov^ing  how 
the  fields  "break  up  into  circularly  polarized  components. 

"Finally,  N  is  given  "by  the  first  line  of  eqviation  4.8),  where  y_ 

«j  J 

appears  in  equation  2.11).  Now  N_,  IL,,  gnd  1T„  may  he  interpreted  as  follows. 
The  incident  field,  velocity,  and  a.c,  charge-  density  are  presumed  to  be  known. 
The  effect  of  incident  velocity  and  charge  is  siammed  up  by  N^,  which  says, 
form  the  linear  combin?tion  F_,  which  gives  the  relative  effect  of  velocity 
and  charge.   Multiply  by  I  (iy  r)/l  (•»t,b),  the  transfer  coefficient  giving 
the  CTirrent  down  the  tube  due  to  a  unit  current  at  the  point  r,  and  integrate 

over  r  .  Then  multiply  by  Jk  { Ij^(p3^b)-G^(p^a)K^(p  b)} /(s^^+Jk),  and  divide  by 

'  /   \  s-i  z 

M,  (s^)  to  obtain  the  actual  current  in  the  mode  whose  z-dependence  is  as  e   , 

Similarly,  IL,  and  N  give  the  effect  of  incident  transverse  magnetic  and  trans- 
verse electric  fields,  respectively,  and  the  appropriate  transfer  coefficients 
show  what  current  is  generated  by  nnit  input  field.  For  values  of  n  other 
than  zero,  the  equation  G  (prO  =  0  may  have  several  roots.  However,  in  gen- 
eral these  will  not  be  in  step  with  the  beam.  Also,  the  contribution  of  these 
higher  modes  is  small  if  the  beam  radius  is  small  compared  to  the  wavelength 
in  the  input  waveguide.  We  now  briefly  consider  the  continuous  spectrum, 
which  arises  from  the  integrals  over  the  branch  cuts.  The  integral  along 
the  upper  eno   lower  sides  of  the  negative  real  axis  may  be  regarded  as  a 
superposition  of  attenuated  waves,  since  s  is  real  and  negative  on  this  con- 
tour. Such  a  superposition  will  generally  have  an  approximate  l/z  dependence 
for  large  7..      It  accordingly  may  be  neglected  when  the  length  of  the  helix  is 
large . 
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The  integral  over  the  part  of  the  imaginary  axis  from  -Jk  to  Jk  is  a 

superposition  of  fast  propagated  waves.  Because  of  the  "boundedness  of  the  inte- 

-1/2 

grand  along  the  contour,  this  integral  will  he  of  no  lower  order  than  z  '   for 

large  z,  and  will  he  rapidly  oscillating  in  z.  The  fast  waves  are  not  amplified 

so  this  integral  also  may  he  neglected. 

Finally,  the  integral  along  the  section  of  the  imaginary acis  from 

-j(k  -h)  to  -J(k  +h)  presents  a  more  complicated  prohlem.  This  is  because  of 

I 
the  involved  nature  of  the  singularity  of  the  integrand  at  s  =  -Jk  .  However, 

this  integral  also  involves  only  propagated  waves,  and,  although  no  detailed  proof 
has  heen  given,  we  shall  presume  that  at  large  distances  down  the  tuhe  this  last 
integral  may  also  he  neglected  compared  to  the  amplified  wave. 

To  summarize,  we  have  developed  a  general  theory,  which  gives  the  fields 
down  the  tube  in  terms  of  the  input  fields,  velocity,  and  charge  density.  The 
major  result  is  equation  ^.7),  which  gives  the  coefficient  of  the  amplified  com- 
ponent of  current.  The  application  of  this  theory  to  particular  types  of  excita- 
tion and  the  comparison  of  these  results  with  earlier  methods  will  he  given  in  a 
second  report. 


Date  Due 


PRlNTECl     IN  U.  S.  A. 


:tyu 

15        Sollfrey. 

c.    1  Effect   of  irJ.tial   ccnditiii 
or.  traveling  wave  tubes. 


15 

Sollfre^^ 

ctof  init 

,ial 

^^^cnax  L,ions   on 
wave   tubes. 

trav 

Bit 

-ng 

DATE    -Ji 

SOHHOAE-^S    ^ 

(AVE 

??OOi< 

N.  Y.  U.  Institute  of 
Mathematical  Sciences 

25  Waverly  Place 
New  York  3,  N.  Y. 


